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The dependence of the Raman spectrum on the excitation energy has been investigated for ABA-and ABC- 
stacked few-layer graphene in order to establish the fingerprint of the stacking order and the number of 
layers, which affect the transport and optical properties of few-layer graphene. Five different excitation 
sources with energies of 1.96, 2.33, 2.41, 2.54 and 2.81 eV were used. The position and the line shape of the 
Raman 2D, G*, N, M, and other combination modes show dependence on the excitation energy as well as the 
stacking order and the thickness. One can unambiguously determine the stacking order and the thickness by 
comparing the 2D band spectra measured with 2 different excitation energies or by carefully comparing 
weaker combination Raman modes such as N, M, or LOLA modes. The criteria for unambiguous 
determination of the stacking order and the number of layers up to 5 layers are established. 



Few-layer graphene (FLG) is attracting much interest recently. The electronic band structure is sensitive to 
interlayer interactions which in turn depend on the number of layers and stacking orders. In natural 
graphite, two stable stacking orders, ABA or Bernal stacking and ABC or rhombohedral stacking, exist 13 . 
Even for the same number of layers, different stacking sequences result in different physical properties because of 
the interlayer interactions. Several studies have elucidated the influence of the number of layers and the stacking 
order on transport 4,5 and optical properties 610 . These results suggest that the number of layers and stacking orders 
are crucial to understand the properties of the FLG system. Therefore, identifying these structural differences is 
very important. Supplementary Fig. SI summarizes the phonon modes in ABA- and ABC-stacked 3-layer 
graphene (3 LG). Similar assignments can be made for thicker FLG. The stacking order affects the phonon 
dispersion through different interlayer interaction. Raman spectroscopy has been established as the most power- 
ful tool to determine the structure of FLG nondestructively. The line shape of the Raman 2D band has been used to 
determine the number of layers 6-8,10,11 and stacking orders 9,12 . Early studies 6 8,10,11 established Raman fingerprints 
to distinguish the number of layers in FLG using the 514.5 or 532-nm lasers without specific reference to the 
stacking order. A later study by Lui et al. showed that Raman imaging technique can be used to visualize the ABA 
and ABC domains on 3 LG and 4 LG samples'. Because of the double-resonance Raman scattering process for the 
2D band, its line shape reflects not only the phonon dispersion but also the electronic band structure, which in 
turn are determined by the stacking order and the number of layers. However, as we show in this work, the line 
shapes of the 2D band from FLG samples with different stacking orders or numbers of layers are sometimes very 
similar for a given excitation energy, making it difficult to determine these structural parameters unambiguously. 
On the other hand, the low frequency shear mode of FLG, the so-called C peak 13 , directly reflects the interlayer 
coupling and so shows a clear correlation with the number of layers. However, because of its proximity to the laser 
line, it is very difficult to observe with conventional confocal Raman apparatuses. Recently, several studies 
demonstrated that other resonant Raman modes, such as so-called N mode 14,15 , M mode 1418 , or other combina- 
tional modes 17,18 also reflect the stacking order and the number of layers. One may, therefore, determine the 
stacking order and the number of layers unambiguously, by combining information from several such modes. It 
should be pointed out, however, that these resonant modes are dispersive because the momenta of the phonons 
involved are dependent on the excitation energy. Therefore, it is essential to have comprehensive understanding 
of the dependences of these Raman features on the excitation energy as well as the stacking order and the number 
of layers. As different researchers use different lasers for the Raman measurements, the excitation energy 
dependence is a crucial piece of information for using the Raman spectrum as the fingerprint of the stacking 
order and the number of layers. 
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In this work, we used five most popular excitation laser wave- 
lengths (632.8, 532, 514.5, 488, and 441.6 nm, corresponding to 
1.96, 2.33, 2.41, 2.54, and 2.81 eV, respectively) to investigate various 
Raman features of FLG. It is found that typically 2 different excitation 
energies are needed to determine the structure by using only the 2D 
band. The G* band is quite useful for differentiating stacking orders 
if a high-energy excitation laser is used. Other weaker modes also 
depend sensitively on the structural parameters if an appropriate 
excitation energy is used. Our work shows how the stacking order 
and the number of layers can be determined unambiguously from 
Raman measurements. 

Results 

2D band. The 2D band, also known as the G' band, originates from 
the double-resonant process which involves two transverse optical 
(TO) phonons near the K points in the Brillouin zone 19 Figure 1 
shows the Raman 2D bands of FLG with ABA- and ABC- stacking 
orders under 5 different excitation energies. The excitation energy 
dependence of the 2D mode without regard to the stacking order has 
been studied by several groups 6,7-20 . Cong et al. 12 and Lui et aV studied 
ABA- and ABC-stacked 3 LG with different excitation energies. In 
our work, we investigated the excitation energy dependence for 1 to 5 
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Figure 1 | 2D band of FLG with 2 different stacking orders under 5 
different excitation energies. The blue, cyan, green, broken-green, and red 
lines represent the measured Raman spectra using 5 different excitation 
energies of 2.81, 2.54, 2.41, 2.33, and 1.96 eV, respectively. The 2D peak 
intensity is normalized to the G peak intensity for each spectrum. With 
increasing excitation energy, the 2D band blueshifts and changes its line 
shape. The line shapes of the 2D band are distinct between ABA- and ABC- 
stacking orders for all excitation energies. 



layer graphene (1 LG-5 LG) with both ABA- and ABC-stacking. 
The blueshift of the 2D bands with increasing excitation energy is 
explained by the resonant scattering process which involves phonons 
with momenta in the TK direction 21,22 . The line shape of the 2D band 
reflects the contributions of scattering paths involving different 
valence and conduction bands, and so the frequency of the 
components reflect the details of the band structure near the given 
excitation energy. Furthermore, the electron-phonon interaction 
could in principle depend on the energy and momenta of the 
electrons and phonons, which would lead to excitation energy 
dependence of the relative intensities of different components of 
the 2D band. 

Due to the dramatic influence of the number of layers and the 
stacking order on the electronic band structure, the line shape of 
the 2D band exhibits clear dependence on the number of layer, the 
stacking order, and the excitation energy and therefore can be used as 
the fingerprint of the number of layers and the stacking order if 
appropriate excitation energy is chosen. For example, with the 
2.54 eV excitation, most cases can be determined with a possible 
exception of ABA-stacked 4 (4 LG) and 5 layer graphene (5 LG). 
In order to distinguish them, one would need to use a lower energy 
excitation source such as 1.96 or 2.33 eV. On the other hand, if one 
uses the 1.96 eV excitation, all but the ABC-stacked 3 and 4 layer 
graphene can be identified. To distinguish them, one needs to use 
higher energy excitations such as 2.54 or 2.81 eV. 

G* band. The G* band which appears in the 2400-2450 cm 1 region 
of the spectra is interpreted as a combination of a transverse optical 
(TO) and longitudinal acoustic (LA) phonon mode near the K point 
resulting from double resonance Raman scattering processes 23 and is 
commonly denoted as D + D". Figure 2 shows the G* band of our 
graphene samples. It has an asymmetric line shape which can be 
fitted with 2 Lorentzian components. It has been suggested that 
the main peak originates from phonons along the KT high 
symmetry line and that the asymmetry is due to phonons with 
momenta away from that high symmetry line 23 . This band shifts to 
the lower frequency as the excitation energy is increased from 1 .96 to 
2.81 eV. 

Figure 2 clearly shows the dependence of the line shape on the 
number of layers and the stacking order. The dependence is stronger 
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Figure 2 | G* band of FLG for excitation energies of (a) 1.96, (b) 2.33, (c) 
2.41, (d) 2.54, and (e) 2.81 eV. The solid spectra are from monolayer and 
ABA-stacked graphene. The dotted spectra are for ABC-stacked FLG. The 
G* band redshifts when excitation energy increases. We observed the 
difference in G* band between ABA- and ABC-stacked FLG when 
441.6 nm or 488 nm lasers were used. 
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Figure 3 | M band in FLG measured under 5 different excitation energies. 
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for higher excitation energies. For example, with the 1.96 eV excita- 
tion, it is difficult to distinguish 1 to 3 layer graphene, and the dif- 
ference between ABA- and ABC stacked FLG is minimal. However, 
with the 2.81 eV excitation, single- and bi-layer graphene can be 
easily identified with respect to thicker FLG. Furthermore, ABA- 
and ABC-stacking can be distinguished for 3 to 5 layer graphene, 
although for the same stacking, the number of layers cannot be 
determined conclusively. In general, the ABC-stacked FLG has smal- 
ler intensity on the lower frequency side of the G* band. The lower 
frequency component is due to phonons along the KT high sym- 
metry line according to Ref. 23. Therefore, its smaller intensity can be 
interpreted as being due to weaker electron-phonon coupling of this 
particular phonon mode in ABC-stacked FLG. 

M band. In the range of 1700 to 1800 cm -1 , a series of Raman peaks, 
called the M band, appear as shown in Fig. 3. The intensity and the 
number of peaks depend strongly on the stacking order and the 
number of layers. These peaks have been interpreted as 
combinations of an LO and an out-of-plane breathing mode ZO' 
(LO + ZO') and an overtone of another out-of-plane breathing 
mode ZO (2ZO) 1418 ; or as overtones of out-of-plane transverse 
optical oTO modes (2oTO's) 1617 . This band is not observed for 
1LG as other researchers reported 141618 . For the ABA-stacked 
FLG, a sharp peak dominates with smaller side peaks. The line 
shape does not change substantially with the number of layers for 
3 or more layers. For the ABC-stacked FLG, on the other hand, there 
are several sharp peaks, and the number of peaks increases with the 
number of layers. So this band can be used to distinguish ABA and 
ABC stackings and also to identify the number of layers for ABC- 
stacked FLG. The positions of these sharp peaks have definitive 
dependence on the excitation energy as shown in Fig. 4. At lower 
energy excitations, the peak separations are larger, and so it is easier 
to distinguish the stacking orders. 

Weak combination modes (1780-2250 cm -1 ). Recently, several 
weak features have been discovered in the Raman spectra of FLG 
in the range between 1780 and 2250 cnT 1 14 ' 16 ls . These features 
originate from combinations of Raman modes which result from 
double-resonance Raman scattering processes involving LO, LA, 
in-plane transverse optical (iTO), and in-plane transverse acoustic 
(iTA) modes. These peaks are shown in Fig. 5. 



For the same type of stacking the spectrum does not vary much 
with the number of layers. On the other hand, there is a clear differ- 
ence between the spectra for ABA-stacked FLG and ABC-stacked 
FLG. The LOLA mode in particular can be used as a fingerprint of the 
stacking. This mode is sharper and appears at higher frequencies in 
ABC-stacked FLG. This is similar to what Cong et al. observed for 
3 LG 12 . The iTOTA (iTO + iTA) mode also shows different line 
shapes for the two stackings, although the difference is rather subtle. 
These trends are observed for all excitation energies used in our study 
(See Supplementary Fig. S7 online). 

The correlation between the peak position of the combination 
modes from 1780 to 2250 cm" 1 and the number of layers in ABA- 
and ABC-stacked graphene under different excitation energies is 
summarized in Fig. 6. As the excitation energy increases from 1.96 
to 2.81 eV, the iTALO (1854-1867 orT 1 ) and LOLA peaks (2011- 
2018 citT 1 ) blueshifts. On the contrary, the iTOTA mode redshifts 
with increasing excitation energy. From the excitation energy 
dependence we can get insight upon the scattering mechanism. For 
intra-valley scattering, the momenta of phonons are close to the T 
point and increase as the excitation energy is increased. When an 
acoustic phonon is involved, the strong dispersion of the acoustic 
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Figure 4 | Dependence of the peak positions in the M band of ABA- and 
ABC-stacked FLG on the excitation energy. 



SCIENTIFIC REPORTS | 4:4630 | DOI: 10.1038/srep04630 



3 




ABA 



ABC 



1800 2000 2200 

Raman Shift (cm 1 ) 

Figure 5 | Combination modes in the range of 1700 to 2250 cm -1 for the 
excitation energy of 2.33 eV. Green spectra are for single-layer and ABA- 
stacked FLG. Red spectra are for ABC-stacked FLG. ABA- and ABC- 
stacked FLGs show different M-band peaks and LOLA-mode peaks. 
Spectra taken with other excitation energies are shown in Supplementary 
Fig. S7. 

phonons near the T point would lead to a blueshift of the Raman 
features with increasing excitation energy. This is what is seen for 
LOLA and iTALO modes. For inter-valley scattering, on the other 
hand, higher excitation energy would lead to smaller phonon 
momenta, moving the momenta of the phonons away from the K 
point. For the iTOTA mode, the acoustic phonon frequency 
decreases, whereas the iTO frequency increases because of the 
Kohn anomaly, as the momenta moves away from the K point. 
Since iTA has a steeper dispersion near the K point than iTO, the 
overall effect should be a redshift, as seen in Fig. 6. 

N band. On the lower frequency side of the G peak (—1580 cm -1 ), a 
weak feature, called the N band is observed. Herziger et al. attributed 
these peaks to a combination mode of a Stokes-scattered LO phonon 
and an anti-Stokes scattered ZO' phonon 14 . The stacking- order 
dependence of these peaks was not investigated. Figure 7 shows the 
N band for ABA- and ABC-stacked FLG for 5 different excitation 
energies. The band redshifts with increasing number of layers. As the 
excitation energy increases the difference between ABA- and ABC- 
stacking becomes more dramatic, whereas the peaks are better 
resolved with lower energy excitations. Low-energy excitations can 
be used to determine the number of layers regardless of the stacking 
order, whereas high-energy excitations should be used to distinguish 
the stacking order. 

Discussion 

Although the 2D band is very sensitive to both the number of layers 
and the stacking order, a single excitation energy is not sufficient to 
provide a definitive fingerprinting for ABA- and ABC stacked FLG 
up to 5 layers in thickness. By combining measurements of the 2D 
band spectra with 2 different excitation sources, a low energy laser 
and another high energy one, we can definitively determine the 
stacking order and the number of layers of FLG. Additionally, a 
distinction between ABA- and ABC-stacking is revealed in the G* 
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Figure 6 | Dependence of the peak positions of the combination modes 
on the number of layers in ABA-stacked graphene in the range of 1780 to 
2250 cm -1 . The blue, cyan, green, dotted- green and red lines represent the 
measured peak positions of combination modes using 5 different 
excitation energies of 2.81, 2.54, 2.41, 2.33 and 1.96 eV, respectively. The 
arrows indicate the shifting direction of combination-mode peaks when 
the excitation energy increases. 

band for high-energy excitations. Weak combination modes such as 
N, M, and LOLA bands also exhibit clear dependence on the struc- 
ture of FLG and can be used to confirm the identifications of the 
stacking order and the number of layers. For example, the LOLA- 
mode peak, which blueshifts with excitation energy in the case of 
ABC-stacked FLG only, is useful for determination of the number of 
layers for ABC-stacked FLG. In conclusion, we have shown that 
Raman spectroscopy, with appropriate choice of excitation energy, 
can be used to determine unambiguously the stacking order and the 
number of layers of FLG up to 5 layers. 

We should note that doping of the graphene samples may affect 
their Raman features. It is well known that the G band has a strong 
dependence on the doping density 24,25 . Therefore, the weak Raman 
features that we use to determine the number of layers and the 
stacking order may be affected by the residual doping of the graphene 
samples. However, we did not observe any appreciable difference in 
these Raman features among several samples that we measured. We 
hence conclude that the uncertainty in the determination of the 
number of layers and the stacking order due to residual doping is 
minimal. 



Methods 

Few-layer graphene samples were prepared on Si substrates covered by a layer of 
100 nm oxide (Si0 2 ) by the mechanical exfoliation method. The number of layers was 
identified by optical microscope images, atomic force microscopy, and by the line 
shape of the Raman 2D band. Stacking orders were identified by comparing the 
Raman spectra with previously reported data in the literature 9,12,18 . Areas with dif- 
ferent stacking were identified by the Raman imaging technique (See Supplementary 
Information). 

For the Raman measurements, we used five different excitation sources: the 441.6- 
nm (2.81 eV) line of a He-Cd laser, the 488-nm (2.54 eV) and 514.5-nm (2.41 eV) 
lines of an Ar ion laser, the 532-nm (2.33 eV) line of a diode-pumped frequency- 
doubled solid state Nd:YAG laser, and the 632.8-nm (1.96 eV) line of a He-Ne laser. 
The laser beam was focused onto the sample by a 50X microscope objective lens (0.8 
N.A.), and the scattered light was collected and collimated by the same objective. For 
the 441. 6-, 488-, 514.5- and 532-nm excitation wavelengths, the scattered signal was 
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Figure 7 | N band in FLG under excitation of 5 different lasers. The broad background spectrum, which comes from the Si/Si0 2 substrate, is subtracted 
in the cases of 1.96, 2.33, 2.41, and 2.81 eV. In the cases of 2.54 eV, the background (black spectrum) was not removed so that the N band is better 
resolved. 



dispersed with a Jobin-Yvon Triax 550 spectrometer (1800 grooves/mm) and 
detected with a liquid-nitrogen-cooled back- illuminated charge -coupled -device 
(CCD) detector. For the 632.8-nm excitation laser wavelength, the scattered signal 
was dispersed with a Jobin-Yvon Triax 320 spectrometer (1200 grooves/mm) and 
detected with a thermoelectrically cooled, back- illuminated, deep -depletion CCD 
detector. For low-intensity peaks such as the N band, long integration times up to 40 
minutes were used in order to get spectra with good signal-to-noise ratios. The laser 
power was kept below 1 mW to avoid local heating and damage to the samples. For 
atomic force microscopy, we use a commercial system (NT-MDT NTEGRA Spectra). 
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